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ABSTRACT 

By analogy with recent arguments concerning the mean velocity profile of wall- 
bounded turbulent shear fiows, we suggest that there may exist corrections to the 
2/3 law of Kolmogorov, which are proportional to (InRe)"^ at large Re. Such cor- 
rections to K41 are the only ones permitted if one insists that the functional form 
of statistical averages at large Re be invariant under a natural redefinition of Re. 
The family of curves of the observed longitudinal structure function Dnir^He) 
for different values of Re is bounded by an envelope. In one generic scenario, close 
to the envelope, Dii{r, Re) is of the form assumed by Kolmogorov, with correc- 
tions of 0((lnRe)"^). In an alternative generic scenario, both the Kolmogorov 
constant Ck and corrections to Kolmogorov 's linear relation for the third order 
structure function Dm{r) are proportional to (InRe)"^. Recent experimental 
data of Praskovsky and Oncley appear to show a definite dependence of Ck on 
Re, which if confirmed, would be consistent with the arguments given here. 
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1. Introduction 



The term fully- developed turbulence traditionally refers to a unique state of 
turbulent behaviour believed to occur for sufficiently large but finite Reynolds 
number (Re). This state is characterised by local isotropy and homogeneity and 
associated universal behaviour of statistical properties, such as moments of the 
longitudinal velocity difference Vr = (v(x + r) — v(x)) -r/lrl. The first theoretical 
description along these lines was given by Kolmogorov and Obukhov in 1941 
(referred to as K41). The assumption that such a limiting state exists, and may 
be found at large but finite Re, is non-trivial, and in our view, has not properly 
been established experimentally. The purpose of this note is to investigate how 
a breakdown of this assumption would be manifested. 

The mathematical expression of the assumption of a limiting state of fully- 
developed turbulence is that statistical averages of the fiow exhibit complete 
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similarity with respect to Re. To explain this statement, let us consider the 
second order longitudinal structure function Dnir) = ((fr)^)- In K41, Dii{r) 
has the form 

DLiir) = CK{lrf'\ (1.1) 

where r lies in the inertial range, and e is the mean rate of energy dissipation 
per unit mass. Kolmogorov's form for Dnir) is based upon both dimensional 
considerations, and assumptions about limiting behaviour. Dimensional analysis 
shows that the form of DLL{r) must be given by 

DLL{r) = (er)2/3F(Re, r/L), (1.2) 

where F(x, y) is a universal function to be determined, L is the external or 
integral scale and r is always considered to lie in the inertial range. Kolmogorov 
assumed that in the limits a; — oo and y ^ 0, the function F{x,y) simply takes 
the constant value Ck- In other words, there is complete similarity with respect 
to the variables Re and r/L. 

3—5 

The existence of the limit of F{x, y) as y — > has been questioned due to 
intermittency — fluctuations of the energy dissipation rate about its mean value 
e. Incomplete similarity in the variable r/L would require the non-existence of 
a finite and non-zero limit of F{x, y) as y ^ 0, and leads in the simplest case to 
the form 

DLL{r) = CKierf/' (1-3) 

where a is the so-called intermittency exponent, believed to be small and non- 
negative. 
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In the present note, we argue that there may be an alternative way, which we 
term asymptotic covariance, in which a lack of complete similarity can occur. We 
try to use physical arguments to constrain the mathematical form that this might 
take. Asymptotic covariance in Re would imply that there is no unique limiting 
state of fully-developed turbulence. Instead, the manner in which statistical 
averages evolve with Re for large Re is governed by a functional form that in the 
simplest case is universal. 

Precisely the same set of arguments can be made for a seemingly different, but 
related problem: the mean velocity profile in a wall-bounded turbulent shear fiow. 
There is a well-known (and probably superficial) analogy^ between the boundary 
layers and universal scaling regimes of both the spatial structure of wall-bounded 
turbulent shear flows and the local structure of turbulence. For example, consider 
turbulent flow in a pipe. The viscous wall layer is analogous to the dissipative 
range in fully-developed turbulence; the velocity proflle (conventionally described 
by the universal von Karman-Prandtl logarithmic law) outside the viscous wall 
layer, but on scales much smaller than the pipe radius is analogous to the inertial 
range (conventionally described by K41); and the non- universal flnite-size effects 
on the flow associated with scales of order the pipe radius are analogous to the 
non-universal behaviour of fully-developed turbulence at the integral scale. This 
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problem was already considered in detail, and it was shown that the existing 
data do not exclude the possibility there of asymptotic covariance in Re. 

This paper is organised as follows. In section 2, we review the analysis of 
the wall-bounded shear flow. In particular, we propose a principle of asymptotic 
covariance, in which we insist that the functional form of statistical averages 
at large Re be invariant under redeflnition of Re. This implies that a form of 
incomplete similarity must occur in terms of the variable In Re. In section 3, 
we consider the analogous arguments for the local structure of fully-developed 
turbulence. We conclude in section 4 with a brief discussion of experimental data 
and some flnal comments. 
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2. Wall-bounded turbulent shear flows 



We consider a wall-bounded shear flow which is statistically steady and ho- 
mogeneous in the longitudinal direction. Its properties vary only in the lateral 
direction, perpendicular to the wall. A classic example, which we shall always 
have in mind, is the flow in a pipe far from the entrance and outlet. 

2.1. Mean velocity profile 

, , 8 9 

Von Karman and Prandtl obtained the law for the variation of the mean 
longitudinal velocity u in an intermediate region of the turbulent shear flow, 
outside a small "viscous" sub-layer near the wall. Within this sub-layer, the 
stress due to molecular momentum transfer is comparable in magnitude with 
that due to turbulent momentum transfer by vortices. The universal, Reynolds 
number independent von Karman-Prandtl logarithmic law has the form 

0=-lnr/ + C', (2.1) 

K 

where 

(f) = u/u^, r] = u^y/v, = (r/p)^/^, (2.2) 

y is the distance from the wall, and p and u are the fluid density and kinematic 
viscosity respectively. The constants n (the von Karman constant) and C are 
universal according to the logic of the derivation. The logarithmic law follows 
from a strong assumption of complete similarity, namely that in the intermediate 
region the contribution of the molecular viscosity and the external length scale 
{e.g., the diameter of the pipe) could be completely neglected. 

It was shown in ref. 7 that this assumption is questionable, and an alternative 
relationship was proposed, corresponding to incomplete similarity: 

= Ct, (2.3) 

where the coeflicients C and a may be expressed in terms of the small parameter 

e = t\- (2.4) 

InRe ^ ' 

by the expansions 

a = \e^O{^\ C=-l- + ^ + 0(e). (2.5) 

It is important that the power a in (2.3) depends upon the Reynolds number 
for two reasons. First, it is well-known that a Re-independent power law form is 
inconsistent with the data. 
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Second, Eyink^^ has argued that a Re-independent power law is inconsistent 
with the rigorous bound on the energy dissipation rate given by Doering and 
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Constantin, (making the plausible assumption that their bound can be carried 
over to the pipe geometry), but the form given in (2.5) is consistent with the 
bound. Eyink's argument is as follows: the average dissipation per unit mass e 
is 

L 

e=^Je'{y)dy, e'{y) = u.'^, (2.6) 



which for the von Karman-Prandtl law (2.1) can be estimated to be 

For the velocity profile (2.3), however, it can be shown that 

e'iv) - (^^y "[2-a(a + l)(a + 2)]('+")/(^+°)Re-2a/(i+a)_ ^2.8) 

If a were Re-independent, the dissipation in the intermediate region would vanish 
sharply as Re ^ 00. However, a more physically plausible alternative is that there 
is only a weak dependence of the dissipation on Re for Re — > 00: 

„ n . constant ^, 

Re ~ constant i.e., a ~ — — - — . (2.9) 

InRe ^ ' 

The corresponding result for the dissipation is that 

which is very similar to the estimate for the von Karman-Prandtl law. 

It seems that the scaling law (2.3) manifests a lack of universality, by virtue of 
its dependence on Re. However, this is not so, in the following sense. Instead of 
the traditional universal straight line in the ^-Inr] plane corresponding to (2.1), 
there is a one-parameter family of curves (2.3) occupying a certain region of the 
plane. This region is nevertheless universal, in the sense that it is bounded by the 
envelope of the family, which is a universal curve. The equation for the envelope 
is obtained by eliminating InRe from (2.3), written in the form = F(ln?7, InRe), 
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and the tangency condition dF/dlnKe — 0, and is found to have the universal 
form 

This envelope is in fact close to (2.1) even for moderate values of Inr] and Re, if 
K ~ 0.4 and C = 5.1. For presently unachievably large values of Re and Inry, the 
universal envelope (2.11) assumes the form 

0=^lnr7 + y, (2.12) 

where the coefficient 2/-\/3e 0.424 is close to the generally accepted value of the 
von Karman constant k, but 5e/2 6.79 is larger than the generally accepted 
value of 5.1 — 5.5. 

Thus, even when there is a lack of similarity, it seems that for In Re ^ 1 in 
the intermediate region of the shear flow, the state is described by the relation 
(f) — F(ln?7, In Re), universal in the sense that the same function F applies to all 
turbulent shear flows. 

2.2. Asymptotic covariance 

It is instructive that the universal relation (2.3) - (2.5) contains Re only 
through the dependence on In Re. In fact, this is inevitable, and a consequence of 
what we propose to call asymptotic covariance. Asymptotic covariance provides 
a general constraint on the way in which lack of complete similarity in Re may be 
exhibited in a turbulent flow. Let us assume that we consider a simple turbulent 
flow, such as that in a pipe, and that a putative state of isotropic, homogeneous 
turbulence is present on some scale r <^ L, where L is taken to be the diameter of 
the pipe. The nature of the turbulent state should be insensitive to small changes 
in the cross-sectional average input and output flow rate U or the diameter L. 
Thus, we require that at sufficiently large Re, the functional form of statistical 
averages, such as Dnir) or 0, should not be influenced by a redeflnition of Re 
such that 

Re ^ Re' = ZRe = Re + 5Re, 5Re/Re ^ 0. (2.13) 

For example, at sufficiently large Re, we could use for L the radius of the pipe 
instead of the diameter, or for U the maximum velocity instead of the average 
one etc., without changing the functional form of Dii{r) or 4>. Let us suppose 
that we consider such a statistical average, whose dependence on Re (and 77 in 
the case of the shear flow considered in this section) is through its dependence 



6 



on a function F(ln?7, -(/'(Re)) having for definiteness a uniformly bounded first 
derivative, and with an unboundedly growing function of its argument, to be 
determined by the following considerations. 

For any Z > 

Re 



V^lRe) = V(Reo) + / ^^^dRe' 
J aRe 



(2.14) 

Re ^ ^ ^ 

V'(^Re) = V(^Reo) + / ^^l^MrfRe'. 

J aRe 

Reo 

Here, Reo is some reference value of Re. The first term on the right hand sides 
of (2.14) at large Re is small in comparison with the second term, because -0 is 
unbounded. However, d'ilj{ZKe)/dKe = Zi\)\Z^<d) where ' denotes differentia- 
tion with respect to the argument. Asymptotic covariance is equivalent to the 
statement that 

^A'(Re) = ^A'(ZRe). (2.15) 
Thus, the right hand side of (2.15) is simply Zip'{ZRe), and 

ReV^'(Re) = ZReV''(^Re) (2.16) 

giving if; (Re) oc In Re. 

Let us now consider the question of whether one could use the Taylor mi- 
croscale Reynolds number ReA, which is often assumed to vary approximately as 
-\/Re. From the point of view of advanced similarity methods and dimensional 
analysis, the important point about Re is that it represents a characterization of 
the system which can be made a priori. That is, it is not an emergent property 
of the flow {i.e., depending upon the solution of the equation of motion), but a 
property of the constraints or boundary conditions placed upon the flow. On the 
other hand, Re^ represents the response of the flow, and in general, for arbitrary 
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Re, will not necessarily have a unique, universal dependence on Re. However, 
for an asymptotically covariant theory, at very large Re, the variation of InRe^ 
with In Re may be written in the form 
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where ao, ai etc. are constants, so that 

ReA ~ Re'*°(lnRe)''i (2.18) 
showing that there is a unique power law relationship between Re^ and Re. 

3. Local structure of fully- developed turbulence 

As explained in the Introduction, Kolmogorov's form (1.1) for 

DlUt) = (er)2/3F(Re, r/L) (3.1) 

is based upon both dimensional considerations, and assumptions about limiting 
behaviour. Let us now examine the question of the lack of complete similarity of 
(1.2), based upon considerations analogous to those used above for wall-bounded 
turbulent shear flows. 

We will assume that the spatial fluctuations in the energy fliix diminish at 
higher Re, because the flow configurations or processes which correspond to these 
fluctuations become increasingly dense throughout the flow. In this picture, the 
scaling assumed in K41 becomes more and more accurate at very high Re. We 
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will refer to this assumption as asymptotic Kolmogorov scaling. 

The behaviour of F{x, y) as x — > oo is the central issue on which we focus. The 
manner in which F could fail to attain a finite limit is constrained by asymptotic 
covariance, which implies that 

DlUt) = (er)2/3G'(lnRe, r/L), (3.2) 

for some function G to be determined, and intermittency, i. e. incomplete similar- 
ity of G{x, y) in the limit y — > 0. Conventionally, it is assumed that the manner 
of violation of K41 is that G{x, y) '-^ y" as y ^ 0, where a is an intermittency 
correction, assumed to be Re-independent and estimated experimentally to be 
positive. In contrast, we wish to investigate the consequences of assuming that 
as a; ^ oo and y — > 

o/Q /7.xQ(lnRe) 
DlUt) = (6r)2/3^(lnRe) (^) (3.3) 

where the prefactor A and intermittency exponent a depend on In Re. In partic- 
ular, the assumption of asymptotic covariance implies that we can write 

a = ao + + 0((lnRe)-2); (3.4) 

In Re 

further assuming asymptotic Kolmogorov scaling implies that ckq = 0. Such a 
form is consistent with the experimental results of Castaing et al.^^ As before. 
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instead of a universal straight line in the Dii{r) / (er)^/"^ - (r/L) plane, we obtain 
a one parameter family of curves, occupying a certain portion of the plane. The 
boundary of this family is a universal, Re-independent curve, which satisfies both 
(3.3) and the the condition dDLL{r)/dRe = 0, i.e. 

dA , ^ ,0:1 In (r/L) , , 

= ^(lnRe) (3-5) 



ci(lnRe) ' ' (InRe)^ 

Empirically it is found that the intermittency correction a is positive, so that 
CKi > 0. However, in the inertial range, r <^ L, so that dA/ diXrvKe) < 0. Thus, 
we may expand A in the small parameter e = 1 /In Re 

A{\nRe) = Ao + Aie + 0{€-^), (3.6) 

where Aq, Ai are non-negative constants. There are two cases to consider: (a) 
Aq 7^ 0, and (b) Aq = 0. We will briefly discuss the experimental estimates of 
Ck in the following section, but for now, we remark that both the scatter and 
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the recent report of a systematic dependence on Re encourage us to consider 
not only the conventional case (a) but also case (b). 

In case (a), both intermittency corrections to the exponent a and to the 
Kolmogorov constant ^0 vanish logarithmically as In Re —> 00. More precisely, 
the relation (3.5) gives 

Qiln jr/L) ^ _^4i 1 

InRe InRe [Ao+^ + 0((lnRe)-2)] ^ " ' 

so that the envelope is given by 

DfTir) = (6r)2/3 (1 + 0(ln(r/L)). (3.8) 

Note the correction term in this formula: the cancellations that occur to first 
order in e make the deviations from K41 very difficult to detect. Indeed, near the 
envelope, individual plots of -Dll(^) at fixed but large Re will exhibit corrections 
to K41 only of order 0{e'^). 

In case (b), the results are more dramatic. The envelope condition (3.5) 
becomes 

"-^ (3.9) 



d(lnRe) InRe 
so that we obtain the condition: 

CKiln (r/L) 



(3.10) 



InRe 

The universal boundary in the Di,L{r)/{er)'^/^ - {r/L) plane is then represented 
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by the curve 

(er)2/3 eailn(r/L)' ^ ' 

The equation of a member of the family of curves at fixed Re, near the envelope 
is 

which is K41, but with a Kolmogorov constant Ck = ^i/elnReA- 

The phenomenological considerations given above also have implications for 
higher moments. For example, consider the third moment DLLL{r) = ((i;^)^), 
which is related to DLL{r) by Kolmogorov's relation 

D,,,ir) = -\,r^,/-^. (3.13) 

Using the form (3.3) for DLL{r), and assuming case (b) above, we find that for 
large In Re 

DLLLir) = -fr (l - ^ (^)'^%-3-/4 + 0((lnRe)-2)) . (3.14) 

The correction to Kolmogorov's "4/5 law" is Reynolds number dependent, in 
contrast to the situation in K41. 



4. Discussion 

Sreenivasan has recently surveyed the available data on Ck, for different 
flows and Re, and came to the conclusion that there is no Reynolds number 
dependence for the data taken as a whole. However, there is considerable scatter 
in the data, and as noted by Sreenivasan, a controlled series of measurements in 
a single flow geometry for a wide range of values of Re has yet to be performed. 

To our knowledge, the only measurements which explore the systematics of 

the possible Reynolds number dependence of Ck are those of Praskovsky and 
17. 

Oncley, in which a weak dependence of Ck on Re is reported. These authors 
attempted to fit this dependence with a power law, but, as shown in figure 1, 
their data are consistent with the logarithmic forms proposed here. Figure 1 
shows the data, together with two best fits to the form given by (3.6), for the 
two cases Aq ^ and = 0. The former case, with a non-zero limit of Ck 
as Re — > oo, may be a slightly better fit than the latter. One measure of this 
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is shown in the insets of figure 1: Inset (A) shows the scatter when C^^-lnRe^, 
is plotted against Re;^, whilst inset (B) shows the scatter when AC^^lnRe;^ is 
plotted against Re^, where AC^^ = Ck — 0.45. A slight upwards trend with 
increasing Re;^ may be discernable in (A), whereas the scatter in (B) seems to be 
more uniform. Clearly, these data (with only eight points) are not suflacient to 
draw any strong conclusions; nevertheless, they are not inconsistent with either 
of the two possibilities ^ or = in (3.6). 

Therefore, it would be of considerable fundamental value to obtain a con- 
firmation (or otherwise) of the results of Praskovsky and Oncley, with greater 
precision. In particular, evidence that = would indicate that fully-developed 
turbulence is not a unique, Reynolds number independent state, approached at 
large enough Reynolds number, but instead is a state with universal Reynolds 
number and external scale dependence, identical for different flows, but with no 
attainable limit. 

In conclusion, we have made two main points in this paper. First, we have 
suggested that there may not be complete similarity in Re — in particular, we 
have suggested that the lack of complete similarity is exhibited as asymptotic 
covariance — and we have detailed some of the consequences of this. In particular, 
a universal scaling should be approached at large Re. Second, the universal 
scaling may be that of K41, indicating the existence of a unique fully-developed 
turbulent state, or may have no asymptotic limit, indicating the non-existence of 
such a simple state. 
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FIGURE CAPTIONS 



1. The Kolmogorov constant, as measured by Praskovsky and Oncley, plot- 
ted on a logarithmic scale against Taylor microscale Reynolds number. Also 
shown are best fits to the functional form given by (3.6), for the two cases 
Aq ^ and Aq = 0. Inset (A) shows the scatter when Ct^IuRca is plot- 
ted against Re;^- Inset (B) shows the scatter when ACi^lnRe;^ is plotted 
against Re a, where ACk = Ck — 0.45. 
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